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ABSTRACT 

Using an analogy between the conductivity tensor of electronic systems and the spin 
stiffness tensor of spin systems, we introduce the concept of the Thouless number gQ and 
the dimensionless frequency dependent conductance g{uj) for quantum spin models. It is 
shown that spin diffusion implies the vanishing of the Drude peak of g{io), and that the 
spin diffusion coefficient Dg is proportional to gQ. We develop a new method based on 
the Thouless number to calculate Dg, and present results for Dg in the nearest- neighbor 
quantum Heisenberg ferromagnet at infinite temperatures for arbitrary dimension d and 
spin S. 

1. Introduction 

According to Thouless0 the conductivity of an electronic system of linear size L is a 
measure of the rigidity of the wave functions with respect to a twist in the boundary 
conditions. This is most clearly seen by writing the real part of the zero- frequency 
limit of the dimensionless conductance in d dimensions, 

in the formBi 

Ec 

go = limReg{uj) = . (1.2) 

Here A is the average spacing between energy levels at the Fermi energy, and Ec 
is the Thouless energy, which is the typical fluctuation in energy levels caused by 
replacing periodic with antiperiodic boundary conditions. In a system with a finite 
diffusion coefficient D, the Thouless energy is given by Ec — hD/L^, and can be 
interpreted as h divided by the time taken by a particle to diffusej-across a box of 
side L. The dimensionless quantity go is called Thouless numbeia. The metallic 
state is defined by go ^ Ij i-e. an electronic system is a metal if an interval of width 
Ec around the Fermi energy contains many energy levels. 
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In the present work we shaU show that g{uj) has a precise dimensionless analog 
g{uj) in quantum spin systems, which can be very useful for a better understand- 
ing of spin diffusion. For simplicity we shall focus here on the spin-S" quantum 
Heisenberg ferromagnet, but it seems that our results can be generalized to Heisen- 
berg antiferromagnets, and models of itinerant magnetism. The hamiltonian of the 
nearest neighbor Heisenberg ferromagnet is given by 

d 

W = -^EESr-Sr+a, , (1.3) 
r /I— 1 

where the r-sum over the N sites of a d-dimensional lattice, and a^, /i = 1, . . . , d, 
are vectors of length a connecting site r with its nearest neighbor in direction /x. 
We restrict ourselves to a hypercubic lattice, where • = a^J^^. J > is the 
exchange coupling, and the Sr are SU{2) spin operators satisfying ~ S{S + 1). 

In their elegant quantum fluids approach to frustrated quantum antiferromag- 
nets Chandra, Coleman, and LarkinQ recently introduced the notion of the dynamic 
spin stiffness tensor i^^-'j^, which is a tensor both in spin space and in real space. 
Here and below the indices /i, — 1, . . . , d refer to the d directions in real space, 
and i,j = x,y, z refer to the three components of the spin operators. A detailed dia- 
cussion of the physical meaning of the spin stiffness tensor has been given in Ref.Q, 
and will not be repeated here. Roughly, fC^^^ measures the energy change induced 
by a space- and time dependent local twist in the direction of the quantization axis 
of the spins. This definition is a generalization of the static spin stiffness p^i which 
corresponds to a time-independent spiral twist of the quantization axis in the limit 
that the wavelength of the spiral becomes infinitely largJl. The existence of close 
analogies between the spin stiffness of classical Heisenberg models and the conduc- 
tance of disordered electrons has first been noticed by ChakravartjU. In Ref.i we 
have derived a spectral representation for K^^^ in quantum Heiseberg ferromagnets, 
and pointed out a formal similarity with the Kubo formula for the conductivity of 
electrons. Similar to the current response-kernel of an electronic system defined in 
the appendix, the spin stiffness tensor has a diamagnetic- and a paramagnetic part, 

K%(k,E) = D% + P;i(k,E) , (1.4) 

where and h/E are the wavelength and time scale characterizing the local twist 
of the spin directions. To make the analogy with the electronic problem manifest, 
we define a mass by setting 



msO? 



J . (1.5) 



The diamagnetic part of the spin stiffness tensor can then be written as 
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and the paramagnetic part P*-J,(k, E) has the spectral representation 



n,rn 



< n|J^(k)|m >< m|J^(-k)|n > 



Em — En — E 

< n\Jl{k)\m >< m| J^(-k)|n > 



Er, 



E„, + E 



(1.7) 



Here En and |n > are exact eigenvalues and eigenstates of Eq jl.3| , p„ are the Boltz- 
mann factors p„ — e~^"/-^/Tr[e~^/-^], and = [J'^, Jji, J^] is the spin current 
operator in direction a^, 



JM(k) 



(1.8) 



The temperature T is measured in units of energy. In d dimensions one can define 3d 
spin currents, corresponding to the three spin projections and the d directions in real 
space. In a temperature regime where the system has long-range spin correlations, 
the uniform spin stiffness Pq is finitecl. In terms of the spin stiffness tensor defined 
above, Pq is given by 



Pa 



-•^ hm 



lim Kl^{\i,huo + iQ"' 



(1.9) 



In case that the system has a spontaneous magnetization, we shall choose a co- 
ordinate system such that the direction defined by the magnetization is the z- 
direction in spin-space. From Eqs.l.E and 1.7 it is clear that the imaginary part of 
iV'fj^(k, ?iw + iO+) vanishes at zero frequency, to that pg is real. At T = Eq.1.9 



reduces to the famihar result pg ~ a'^^'^JS^. For temperatures above the ordering 
temperature Tc, the diamagnetic- and paramagnetic contributions to Pq precisely 
cancei^'@, so that = 0. In fact, pg is the magnetic analog of the the long- 
wavelength limit of the London-Kernel in an electronic system. Both quantities 
describe the appearance of a long-range rigidity in the wave functions. An obvious 
question, which apparently has not been discussed in the literature, is whether in a 
localized spin model one can also define the analog of the weight of the Drude peak 
in the expression for the conductivity of a metal. Furthermore, what is the analog 
of the Thouless number in spin systems? In the present paper we shall answer these 
questions. 

2. Thouless number and spin conductance 

Guided by the rescaled form of the Kubo formula for electrons given in the appendix, 
we define for the spin system the functions (?(a;), Pii^), Kii-^), and if (k) by simply 

and A. 9 the rescaled current response kernel Ku by 



replacing in Eqs.A.f; ,A.6, A 



the spin stiffness tensor ifff , 
5(w) 



lim 



i{TiLj + iO+) 



(2.1) 
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P(uj) = lim 



ImXff(k, huj + iO^ 

huj 



lim Re/Cff (k, huj + iO^ 



K'ik) = - lim K^^{k, huj + iO+) 



(2.2) 
(2.3) 
(2.4) 



We also define spin analogs Kq and Kq of the weight of the Drude-peak and the 
long-wavelength limit of the London-kernel, 



A'o = Hm K(uj) 



n = limii:''(k) 
k^o 



(2.5) 
(2.6) 



The function g{uj) plays the role of the dimensionless conductance (see Eq. |A.5| ), and 
we shall call g{LLj) "spin conductance" . P(w) is the paramagnetic contribution to the 
real part of g(w) (see Eq. A.8), and the zero-frequency limit of K{u!) yields the spin- 
analog of the weight of the Drude peak (see Eq. |A.S| ) . The wave- vector dependent 
function K^Ck) is related to the length-scale dependent spin stiffnessQ p'^i'k) = 
L^^'^K'^Ck), and corresponds to the London- kernel of an electronic system (see 



Eq. A.6). A Schwinger-Boson calculation of p*(k) for the two-dimensional quantum 
ferromagnet has been given in RefE. Because the imaginary part of Kfi{k,iO^) 
vanishes, it is not necessary to take the real part on the right-hand side of Eq. ^.4| 
We now connect the above definitions with the phenomenon of spin diffusion. It 
is easy to obtain the following exact spectral representations, 



P(^) 



TT I I lini X!^"*^ ~ - 'hijj)\< n| (k)|TO > 



lim Y''P\- 

k^o [ _ 



Pn Pm 



Em ~ En — hoj 



< n\J^{k)\m > 



(2.7) 



(2.8) 



where V denotes the Cauchy principal value. General hydrodynamic argumentsu 
tell us that spin diffusion can only exist in the paramagnetic regime T > Tc- In 
this case the low-frequency and long-wavelength behavior of the dynamic structure 
factor S{k, uj) is of the form 



S'(k, w) = 2TThJ2PnKEm - E,, - hLo)\ < n\S^\m > 



2X 



(2.9) 



where Ds is the spin diffusion coefficient, and x = T ^ < S^S^ > is the 
uniform susceptibility. The Fourier transform of the spin operators is defined by 
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Sk = X^r G^^ '^Sr- The Heisenberg equation of motion for Sk yields 

d 



q M=i 



(2.10) 

where the momentum sum is over the first Brihouin zone. On the other hand, 
Fourier transformation of the right-hand side of Eq.1.8 yields 



. J a 



Jp(k) = i— Z Sq X Sk-qSin [(k - q) • 



(2.11) 



where we have used ?i^/(TOsa^) = J. Combining Eqs. 2.1C and 2.11 and taking 
matrix elements we obtain to leading order in |k • a^| <C 1 

L 

[E^ ~ E„,) < n|Sk|m >= V < n|J^(k)|m > k • . (2.12) 



Inserting this into Eq.2.7 we arrive at 



P{uj) = N- 



2h 



2 _ ^~hu}/T 



flLO 



5(k, 



lim , 



(2.13) 



we 



Finally, using the assumption of the diffusive form of S'(k, cj) given in Eq.2.9 
obtain 

N 

\imP{uj)^—X^D, , (2.14) 

To complete the analogy with the electronic system, we should also proof that spin 
diffusion implies the vanishing of the weight of the Drude peak Kq defined in Eq j2.5| . 
To show this, we use the following trick: We know that for T > Tc the uniform 
spin stiffness Pq vanishes @3, so that Kq — — limk-^o linicj^o Kfi(k, uj) — 0. Hence 
Kq = Kq — Kq for T > Tc- But the diamagnetic part D^^^ of the spin stiffness tensor 
is independent of momenta and frequency, and does not appear in the difference 
Kq — Kq, so that 



Ko - H 



lim lim — lim lim 



P7i - Pn 



En 



En — huj 



< n| (k)|m > 1^ 



(2.15) 

The contribution from all non-degenerate states En ^ Em cancels on the right- 
hand side of this identity, because in this case the zero-frequency limit is harmless 
and we may interchange the order of the limits. Degenerate states contribute only 
if the limit w — > is taken before the limit k — + 0, because for finite u; we have 
limB„^_E„(p„ -Pn)/{Em ~ En- fluj) = 0. Hcncc 



(2.16) 
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Using the identity 



dw' ^ 5{E,n ~ En ~ huj') 



(2.17) 



E„=E„ 



and comparing the right-hand side of Eq. ^.l6 with the spectral representation of 
P{uj) in Eq. 2.7, we conclude that 



Ko - hm / duj'P{uj') + K^Q 



(2.18) 



Using now the fact that according to Eq. 2.14 P{lu) has a finite limit as u; ^ 0, 
and that — for T > Tc, we see that the right-hand side of Eq.2.18 vanishes, 
so that Kq = 0. The vanishing of the Drude peak in the presence of diffusion is 



familiar from the electronic problem, see Eq.A.12. The Thouless number of a spin 
system with spin diffusion is then 



Ec 

go = Hm Reg{uj) = — ^ 

u—tQ A 



(2.19) 



where A — {Nx)~^ plays the role of the level spacing A of the electronic system, 
and the Thouless energy for the spin system is again defined by Ec = hDg/L^- 
Note that the level spacing at the Fermi energy in an electronic system can be 
written as A = (Nk)^^, where k is the compressibility. Comparison with Eq.1.2 
shows that the only difference between the Thouless numbers in electronic- and 
spin systems is that the compressibility of the electronic system is replaced by the 
magnetic susceptibility of the spin system. 



3. Spin Diffusion 



3.1. General remarks 

First of all, it should be emphasized that there exists no proof of spin diffusion in 
Heisenberg ferromagnets. Our calculations in Sec. are based on the assumption that 
the long-wavelength and low-frequency behavior of the dynamic structure factor is 
of the diffusive form given in EqjlJ. The fact that even for temperatures large 
compared with J the spin diffusion problem is highly non-trivial is most clearly 



seen by writing the spin diffusion coefficient in the form (see Eqs.2.7 and 2.14) 



_7 /-oo Ti-N {e-f^^iSin + e- H)iS{e- H)] 
hDsQ-^ = — lim lim / de r^j ^ , (3.1) 



where P = J/T, and Trjv denotes the trace over the Hilbert space of the A''-site 
Heisenberg model, and the dimensionless operators H = + + H'^ and I are 
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^iven by 

H' = -EE^r^r+a, , ^ = X,y,Z (3.2) 
^ = ;^E^r [^r+a.-^r-aj ■ (3.3) 



d 



From Eq.3.1 it is obvious that a finite spin diffusion coefficient can only be obtained 



in an infinite system, fn any finite system the trace m Eq.3.f will consist of a sum 
of S functions, so that Ds cannot be defined. See Ref.EJ for a discussion of this point 
for the electronic problem. Moreover, even after taking the thermodynamic limit 
— > oo, it is not clear that the right-hand side of Eq jsj] reduces to a finite constant. 
In principle, there are three possibilities, familiar from disordered electronic systems: 



(f) Perfect conductor. Eq.3.1 contains two (5-functions, but only one energy in- 



tegration. If the current operators would commute with the Hamiltonian, then the 



right-hand side of Eq.|3.1| would be proportional to S{^), so that the spin conduc- 
tance has a Drude-peak, and spin diffusion does not occur. In this case the system 
behaves like a perfect conductor. Although in the Heisenberg model [/, H] ^ 0, 
there exists no proof that the non-commutativity is sufficient to remove the Drude 
peak. 

(2) Metal. In this regime the Thouless number go and the spin diffusion coeffi- 



cient are finite, and are related via Eq.2.19 



(3) Insulator. The third possibility is that in a certain parameter regime the 



right-hand side of Eq.3.1 scales to zero in the thermodynamic limit. This would 



correspond to the insulating state of a disordered electronic system. 
3.2. Spin diffusion at infinite temperature 

We now assume that the spin diffusion coefficient is finite and develop a transpar- 
ent and direct method to calculate at infinite temperatures. The problemjif 
calculating at T = oo has been studied intensely more then 20 years ago0-t2l. 
Most methods are based on an indirect calculation of Dg via the dynamic struc- 
ture factor S'(k, w), assuming that its long-wavelength and low-frequency behavior 



is of the form given in Eq,2.£. The spin diffusion coefficient is obtained indirectly 
from S(k, uj) by means of the limiting procedure lim^j^o linrk^o(t^^/fc^)S'(k, us), see 
Eq.^ 

In praxis, it is impossible to calculate S'(k, uj) at low-frequencies, or equivalently 
its real-time Fourier transform ^(k, t) at long times t. The moment method first 
applied by de GenneaJ is equivalent to an extrapolation of a short time expansion 
to long times. The concept of the Thouless number in spin systems offers a more 
direct way to calculate Dg. Of course, if S'(k,aj) could be calculated exactly, then 
the result for would be identical with the result obtained by means an exact 
calculation of g{Lo). However, extrapolations of high-frequency expansions of g{uj) 
and S'(k, u) will in general not agree, because only in the limit w — > and k ^ 
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there is a direct connection between these two quantities. We beheve that our 
approach via the Thouless number 170 is more rehable than extrapolations based 
on the dynamic structure factor, because go is directly proportional to Ds, and no 
further limiting procedures are required. 

We now use the first two expansion coefficients in the short time expansion 
of the right-hand side of Eq. to estimate Ds at T ^ 00. Introducing Fourier 
representations of the (5-functions and defining I{t) — e'^^'/e^*^*, we obtain from 



Eq.3.1 after straightforward manipulations 



7 

msa-' = / C{t) 

Tx Jo 

c{t) = ]-<m\i{t) + i{-t) 



> 



(3.4) 
(3.5) 



where < • • • > denotes thermal average with the Hamiltonian H. The existence 
of the integral in Eq.3.5 implies spin diffusion. The convergence of the integral 
is determined by the long-time behavior of C'{t). Unfortunately, there exists no 
completely controlled method to calculate C{t) for large t. We therefore assume 
that the integral exists, and try to extract the long-time behavior from the short- 
time expansion, 

(-l)"t2n 



(2n)! 



-Co 



(3.6) 



Note that only even powers of t appear, because C{t) = C{—t). The expansion 
coefficients can be written in terms of multiple commutators. The first two 
coefficients are 



Co = <I'> 



C'2 



< I I,H ,H > 



(3.7) 
(3.8) 



where all operators are at equal times. At T = 00 the evaluation of the thermal 
averages in Eqs. [3.7| and |3.8| simplify, because spins at different sites are not corre- 
lated. To calculate Co, we need < {S^Y >= S{S + l)/3, for i ^ x,y,z. A short 
calculation gives 

'S{S + 1) 



Co = 2 



3 



(3.9) 



The evaluation of C2 is tedious but not difficult, ft involves the expectation values 
of up to four spins, which have been tabulated in Ref.EJ. The following averages 
are needed 



S{S + 1) 
15 



[35'(S' + 1) - 1] 



(3.10) 
(3.11) 
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where < . . . >= Tr[. . ■]/{2S + 1), and the trace is over the 25* + 1 states of the 
spin-5 Hilbert space, e*-'*^ is the antisymmetric Levi-Civita tensor. For i j we 
have 



S{S + 1) 
15 

S{S + 1) 
15 



3(3 + 1) 



[5(5 + 1) -2] 



Writing 



-"2 



i.j—x,y,z i.j—x,y,z 

we obtain for a d-dimensional hypercubic lattice 



> 



"5(5 + 1)' 


2 


r/ _n 45(5 + l)_3l 


3 




V 5/ 3 5 



"5(5 + 


1)1 


2 r 


3 






"5(5 + 


!)■ 


2 r 


3 





3\ 45(5+1) 



^2 — ^2 — 



and Cl^ = Cf = Cf ^ = = 0. Hence 



45(5+1) 1 
15 ^ 5 



C. = 



5(5 + 1) 



(4. -2)1^(^^-1 



(3.12) 
(3.13) 

(3.14) 

(3.15) 
(3.16) 
(3.17) 

(3.18) 



We are now ready to extrapolate C(t) to long times. This extrapolation is of course 
not unique. A widely used extrapolation scheme, which we shall follow here, is to 
assume that the first two coefficients are consistent with a GaussianEJ. This leads 
to 

' O2E 



C{t) K. Co exp 



2Co 



(3.19) 



Using the fact that Tx — 5(5 + l)/3 at T = cx3, we finally obtain from Eqs.3^, 3l9| , 
3.18 and 3.19 for the spin diffusion coefficient at T = 00 





'tt3{3 + 1)' 


1/2 r 


a? J 


3 





4d~2 



-1/2 



(3.20) 



45(5+1) 

This is the main result of this section. Note that the term proportional to [5(5 + 
1)]^^ in the second factor of Eq. 3.2C can be viewed as a quantum correction, that 
becomes irrelevant in the limit of large 5. For 5 = 1/2 Eq. 3.20| reduces to 



a? J 



4d- 3 



1/2 ( 0.886 ford=l 
0.396 for d = 2 
0.295 for d = 3 



(3.21) 
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The classical limit of the Heisenberg model should be taken by letting S ~+ oo while 
keeping Jd = JS^ constant. In this limit Dg vanishes. The leading coefficient for 



large S can be read off from Eq. |3.20 
hD, 1 



4(i- 2 



1/2 ^ ( 0.724 for d = 1 

X I 0.418 for d = 2 (3.22) 
0.324 for d = 3 



S 



Eq.3.21 agrees with the results listed in the first column of table V of RefJll. Hence, 
at least to second order in the short time expansion, our method is equivalent 
with Morita's memory function formalisnJlj. However, it seems that our Thouless- 
numbcr approach is physically more transparent, because it does not involve the 
construction of momentum dependent auxiliary quantities. —The results in d = 
3 are also in excellent agreement with Bennett and MartinlHl, who obtained by 
means of an indirect extrapolation method based on the dynamic structure factor 
hDs/{a^J) ~ 0.33^/S(sTT). This gives nD,/{a^J) = 0.29 for S = 1/2, and 
hDs/ia^Jci) = 0.33/5 for 5 ^ oo. 

4. Conclusions and open problems 

In this paper we have used an analogy between the Kubo formula for disordered 
electrons and the dynamic spin stiffness tensor of localized quantum spin models to 
define several new quantities that characterize the dynamic behavior of Heisenberg 
ferromagnets in the paramagnetic regime. The key observation is that, after proper 
rescaling, the dynamic spin stiffness tensor of the spin system and the current 
response kernel of the electronic system can be written in a formally identical way. 
The proper definition of the Thouless number and the dimensionless conductance 
in the spin system directly follow from this analogy. As a first application, we 
have used the Thouless number to develop a simple extrapolation scheme for the 
calculation of the spin diffusion coefficient at infinite temperatures. 

This work opens a number of new directions for further research: A generaliza- 
tion of the concepts developed here to antiferromagnets or Hubbard models seems 
possible, although some technical difficulties might be encountered. The analogy 
between electrons in the presence of disorder and spin models in the paramagnetic 
regime might serve as a useful guide to understand spin diffusion in two dimensions. 
A Schwinger-Boson calculation of Dg in two-dimensional Heisenberg models at low 
temperatures has recently been given by Chubukovt^. His result for Dg in ferromag- 
nets can be easily reproduced if we translate the Drude formula for the dimensionless 
conductance of an electronic system into spin language. By comparing the ex- 
pressions for the diamagnetic tensors -D^^j^ and D^^^ given in Eqs.1.6 and A. 2, we 
conclude that the energy Ti^ N / {mf.L'^) corresponds to h'^NS^/{msL'^) in the spin 
model. Here we have used the fact that for T <C JS^ the spin-correlation length is 
exponentially large compared with the lattice spacingiH, so that the summation in 
Eq.1.6 yields NS^U. Hence, to obtain the Drude result for Ds, we should replace 



h^N/{meL^) n^NS^/{mgL^) = JS"^ in Eq. |A.13| . This gives for the Thouless 
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number in the Drude approximation 



T 



go^JS'- . (4.1) 

Here r is the characteristic hfetime of single-particle excitations at wavelengths 
large compared with the spin correlation length. To determine r, a microscopic 



calculation is necessaryE£l. Combining Eqs.4.1 and 2.19, we conclude that the spin 



diffusion coefficient in d = 2 is at low temperatures given by 



nD,a-' = :^ , (4.2) 



Eq.4.2 agrees exactly with the result of ChubukovEEl (see his Eqs.2.30 and 2.32), who 
performed a diagrammatic calculation within the Schwinger-Boson formalism. The 
factor of JS^ in the enumerator is interpreted by Chubukov as transport coefffcient. 
Note, however, that in = 2 the Thouless number for non-interacting electrons in 
the presence of disorder scales to zero in the thermodynamic limilfi, so that the 
system is an insulator. Assuming that our analogy can also be applied to this case, 
we speculate that for any T > the two-dimensional Heisenberg model corresponds 
to an insulator, so that the spin diffusion coefficient vanishes. In any case, we 
expect that in two dimensions the interaction between the diffusion modes will 



be very important, so that the Drude result given in Eq.4.2 cannot be trusted. 
Non-perturbative methods are necessary to understand spin diffusion in the low 
temperature regime of two dimensional Heisenberg models. 
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Appendix 

In this appendix we shall rewrite the standard Kubo formulaEl for the conductivity 
of an electronic system in a rescaled form. By comparing the rescaled form of the 
current-response Kernel Kfj^^, with the spin stiffness tensor if^^^ defined in Sec, the 
proper definition of the Thouless number and the dimensionless spin conductance 
become obvious. Except for a tilde, we shall use the same symbols as in the spin 
problem to emphasize the close similarity between the rescaled current-response 
kernel and spin stiffness tensor. 

Consider a d-dimensional box of volume L"^ containing N electrons with mass 
TOg- Within linear response theory, the change in the current density Sjfj_(r,t) due 
to a change in the vector potential SA^{r',t') is given by 



Sj,{r,t)^iaL'-yh) 



dt' 



dv'Y,k^,{v,v',t-t')5A,{v',t') 



(A.l) 



where a = e'^/{hc) « 1/137 is the fine-structure constant. For reasons obvious 
below, we have not included the factor aL'^^'^/h into the definition of K^. The 
Kubo formula for the Fourier transform of the current response kerneO yields 



given by 



D^ll' + Pfii^iM-i E), where the diamagnetic- and paramagnetic parts are 



D,. 



< n| J^(k)|m >< to| J^(-k)|n > 



E„ 



En — E 



< n|X(k)|m >< to| J^(-k)|n > 

Em. — En + E 



and the current operators J/i(k) are defined by 



nieL 



dre 



ik-r 



2i 



h.c. 



(A.2) 



(A.3) 



(A.4) 



Here ip{r) is the usual second quantized field operator for the electrons. Note that 
with this rescaling the has units of energy, just like the spin currents defined in 



Eq.1.8. The dimensionless conductance 5 (w) — L"^ ^cr(a;)/(e^/?i), and the (rescaled) 
London kernel K''{'k) are given by 

Kii{k, hLO + iO+) 



lim _ , , 

k^o i{nuo + iQ+) 



(A.5) 



Thouless number and spin diffusion in quantum Heisenberg ferromagnets 13 



(A.6) 



From Eq. A.5 it is clear that the factor aL'^~'^/fi in Eq. A.l is very natural. Note that 
the "w-limit" (where the limit w ^ is taken after the limit fc 0) and the "fc- 
limit" (where the limit fc — s- is taken after the limit a; — > 0) describe very different 
physical properties of the system. The energy Kq = limk-+o (k) is proportional 
to the density of superconducting electrons. Comparing of Eqs.l.E and Eqs. [A.6| , it 
is evident that L^^'^Kq corresponds to in the spin system. Both quantities are 



only finite in the presence of off-diagonal long-range order. From Eq, A.5 we obtain 



Re^(ti;) = ■jTKoS{fiuj) + P{uj) , l-mg{uj) — 



Tiuj 



with 



Piuj) = lim 



ImJ^ii(k, hui + iQ+) 



limReii'ii(k,?iw-f iO+) 

k— »0 



(A.7) 

(A.8) 
(A.9) 



and Kq = lim^^^o K{lo). It is instuctive to evaluate the above quantities in the sim- 
plest possible approximation, where all scattering processes are taken into account 
by introducing a phenomenological lifetime r in the electronic Greens functions. 
At temperatuxes small compared with the Fermi energy, this yields the well known 
Drude resultsEj 



P{u;) 



hN 



hN 



ILOT 
T 



,L2 1 + {uJTf 



1 - 



1 



1 + (cjr)2 



The Thouless number go defined in Eq. L2 is in this approximation given by 

hNr 

90 



(A.IO) 
(A.ll) 

(A.12) 
(A.13) 



Combining Eqs.1.2 and A.13 , and using the fact that the level spacing at the Fermi 
energy in d dimensions can be expressed in terms of the Fermi velocity vp as A = 
'meVp/{dN), we obtain the classical result for the diffusion coefficient D = Vpr/d. 



